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Abstract 

We present modeling of an incompressible viscous flow through a fracture adjacent to a 
porous medium. We consider a fast stationary flow, predominantly tangential to the porous 
medium. Slow flow in such setting can be described by the Beavers- Joseph-Saffman slip. For 
fast flows, a nonlinear filtration law in the porous medium and a non- linear interface law are 
expected. In this paper we rigorously derive a quadratic effective slip interface law which holds 
for a range of Reynolds numbers and fracture widths. The porous medium flow is described by 
the Darcys law. The result shows that the interface slip law can be nonlinear, independently of 
the regime for the bulk flow. Since most of the interface and boundary slip laws are obtained 
via upscaling of complex systems, the result indicates that studying the inviscid limits for the 
Navier-Stokes equations with linear slip law at the boundary should be rethought. 
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1 Introduction 



Coupling between a fast viscous incompressible fracture flow and an adjacent filtration through 
porous medium occurs in a wide range of industrial processes and natural phenomena. The classical 
approach is to model the fracture flow using the lubrication approximation and to replace it by 
an interface condition. Subsequently, it is coupled with a porous medium flow, described for small 
Reynolds numbers by the Darcy's law and by the Forchheimer's law in the case of large Reynolds' 
number. 

Study of the coupling between slow viscous incompressible fracture flow and a porous medium 
was undertaken in [3] and [4]. For the critical fracture width, the interface condition linked to the 
Reynolds' equation from lubrication was found. 

To describe a contact between a porous medium and a large fracture with the width significantly 
larger than the pore size, the following effective slip interface law was established in the seminal 
work by Beavers and Joseph [2], 



where asj is a dimensionlcss parameter depending on the geometrical structure of the porous 
mediumand K is the scalar permeability. v T is the tangential velocity and n is the unit normal 
exterior to the fluid region. Note that in the original version of the law |T]), v T was replaced by the 
difference between v T and the tangential Darcy velocity at the interface. In |18j . Saffman remarked 
that the tangential Darcy velocity at the interface is of order O(K). Then, the slip law without the 
tangential Darcy velocity at the interface (JTJ) became generally accepted. 

The rigorous derivation of the law by Beavers and Joseph through a homogenization limit and 
by constructing the interface boundary layer was done by Jager and colleges in [TU], [TT] and [12] . 
The pressure jump at the interface was studied analytically in [16| and using numerical simulations 
in [BJ. For the review of the results we refer to [13], Q2] and [?]• 

Sahraoui and Kaviany investigated in [19] a flow at the interface between a fracture and a 
porous medium by direct numerical simulations. The interest of this work was in the interface 
laws in presence of large Reynolds' numbers. The interface slip behavior in that case turned out 
to be complex. It was concluded that the flow inertia effects appear independently from the bulk 
nonlinear filtration in the porous medium. If e is a characteristic nondimensional pore size, then for 
longitudinal Reynolds' numbers of order 0(1/ e), numerical simulations indicate that the slip law 
ceases to be linear. The inertia forces at the interface become significant for Reynolds' numbers 
of order O(0. 1/e). Then, the slip coefficient cxbj increases. For the bulk porous medium flow, 
the nonlinear effects become visible only for Reynolds' numbers greater than 0(3/ e). Those 
observations led to a conclusion that asj depends on the Reynolds' number, [14) and [5]. Similar 
conclusion is in |15) . 

However, it seems that a linear slip law, even with the slip coefficient depending on Reynolds' 
number, is not enough for an accurate approximation and that a nonlinear slip law should be 
derived. We will justify it by constructing rigorously an accurate approximation to the velocity field 
and showing that it leads to a quadratic slip law. 

In the present paper we aim to identify a setting corresponding to a nonlinear slip law. We show 
that for a range of values of Reynolds' number and fracture width, the homogenization leads to a 
nonlinear interface law, even though the bulk filtration remains of the Darcy type. To streamline 
the presentation, we focus on a mathematical model in a simple setting. We consider a constant 
driving force, present only in the fracture and, for simplicity, impose periodic longitudinal boundary 
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conditions for the velocity and for the pressure. Such simplification allows to avoid handling the 
pressure field and the outer boundary layers. The general case of nonstationary flows with physical 
boundary conditions and forcing terms will be considered in forthcoming papers. 

The paper is organized as follows: In section [2J we define the problem as a stationary incom- 
pressible Navier-Stokcs flow with Reynolds' number of the order e -7 and the fracture width of the 
order e s . Assuming a relation between 7 and 5, allows us to obtain an approximation which satisfies 
a nonlinear slip law while keeping a linear filtration in a porous medium. In section [3] we 

construct the approximation and prove that it provides a higher order approximation to the original 
problem. 

2 Main result 

2.1 Geometry 

We consider a two dimensional periodic porous medium Q2 = (0, 1) x (—1,0) with a periodic ar- 
rangement of the pores. The formal description goes along the following lines: 

First, we define the geometrical structure inside the unit cell Y = (0, l) 2 . Let Y s (the solid part) 
be a closed strictly included subset of Y, and Yp = Y\Y S (the fluid part). Then, we introduce a 
periodic repetition of Y s all over R 2 and set Yjf = Y s + k, k G Z 2 . Obviously, the resulting set 
E s = U feGZ 2 Y s k is a closed subset of R 2 and E F = R 2 \E S in an open set in R 2 . We suppose that Y s 
has a smooth boundary. Consequently, Ef is connected and E s is not. Finally, we notice that is 
covered with a regular mesh of size e, each cell being a cube Yf, with 1 < i < N(e) = |f22|e _2 [l+o(l)]. 
Each cube Yf is homeomorphic to Y, by linear homeomorphism Ilf , being composed of translation 
and a homothety of ratio 1/e. 

We define Y§. = (n?) -1 ^) and F|, = (nf)" 1 ^). For sufficiently small e > 0, we consider a 
set T e = {k e % 2 \Y§ k C H 2 } and define 

o e = |J Y§ k , s £ = do e , ni = n 2 \o £ = 2 n sE F . 

Obviously, dfl^ = 9^2 U S £ . The domains O e and fi| represent the solid and the fluid part of the 
porous medium f2, respectively. For simplicity, we assume 1/e G N. 

Let < S < 1. We set E = (0,1) x {0}, n{' S = (0,1) x (O,^) and = (0,1) x (~l,s 5 ). 
Furthermore, let e = 0| U E U s /. 

In such geometry, homogenization of the Stokes equation with no-slip boundary conditions on 
S £ leads to Darcy law (see [T], [8], [20] and [21] ). In the presence of inertia, a nonlinear corrections 
to Darcy law arise, as studied in [5]. 



3 



2.2 Position of the problem and the nonlinear slip law 

Let < 7 < 3/2 and let F be a constant. In e we study the following stationary Navier-Stokes 
equation 

-e 7 Av e + (v £ V)v e + Vp £ =fe 1 l { , 2>0} in n £ (2) 

divv e =0 in il e , [ p e dx = 0, (3) 

J01' 5 



v £ = on dQ £ \ \{xi = 0} U {x\ = 1} J , {v £ ,p e } is 1 - periodic in x\. (4) 

Remark 1. We skip here a discussion of modeling aspects. We only mention that e 7 stands for the 
inverse of Reynolds ' number and that the small fracture width e s prevents creation of the Prandtl 's 
boundary layer. 

In order to simplify calculations we take a constant F. It corresponds to a pressure drop. 
Additionally, we assume it only in the fracture £ { S . Let 



W e = {ze H\n E ) 2 , z = on dn e \ ^{zi = 0} U {X! = 1}J and z is 1 - periodic in x ± }. (5) 

The variational form of problem (HJ-Q reads: 

Find v £ e W € , div v e = in £ and p e e L 2 (n e ) such that 

/ e 7 Vv £ Vv3dx+ / (v £ • V)v e ipdx- / p £ div (pdx= Ftpidx, M^p e W £ . (6) 
Jiv Jn? JQe J01- 6 

Theory of the stationary Navier-Stokes equations with homogeneous boundary conditions results 
in existence of the least one smooth velocity field v e g W e , div v e = in f2 e , which solves ([6]) for 
every ip S W 6 , div tp = in Q E . The construction of the pressure field goes through Dc Rham's 
theorem. For more details we refer to the classical Temam's book [2"2"] . 

Now we make assumptions on the parameters 5 and 7. 
(HI) 2 7 < 35, 

(H2) < S < 1 and < 7 < 3/2, 
(H3) 4(5 < 2 7 + 1. 

Now, we can formulate the main result 

Theorem 2. Let us suppose the hypothesis (TIlj-fHSj and letlA 2,e be defined by 



where the boundary layer functions f3 bl and /3 1, are defined, respectively, by J^i[ )- f77P o,nd 
H66]). The constant C\ l < is the stabilization constant for f3^ when yi — > +00. Similarly C\\ is the 



stabilization constant for when 1/2 
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Then, the following estimate holds 

£||VW 2 < e || i2(0 ^ + ||W 2 < £ || LW +£ 1/2 ||W 2 "IUw +e 1 - 5 \\U^\\ L2{&l , Sy < Ce 7 ' 2 - S ~\ (8) 

Remark 3. The rigorous result from Theorem @ showing that U 2 ' £ is of order 0(e 3 ~ s ^ 1 ) on 
£, allows justifying a nonlinear interface law. Contrary to the classical situation, when Saffman's 
modification of the linear slip law by Beavers and Joseph (see JBjj and \lSj) is used, the nonlinear 
interface laws are rarely derived in the literature. However, they are supposed to be appropriate for 
fast flows. 

Setting (5=1 — lr\jY2 and 7 = 3/2 — 77, where < n < 3/2, which fulfills hypotheses (Hi )-(R3 ), 
we obtain on the interface £ 



Vl (s)\ s = -f ^"^(i - cf^-'^c^b - (f ^^"^^'"(f )| 5 



2 e 2 e 



and for the average over the pore face on £ 



Next, for the shear stress we have 



< v 1 (e)\x >= vf f = -|vi^ /12 (l - Cf £ 7 "/ 12 )C? - (i-) 2 ^ 11 "/ 6 < # W (^)| E > ■ (9) 



After averaging over £ wz£/i respect to y\, we obtain 

< _^£l| s <= = ^ £ -l/2+5Wl2(l + e W C U _ e 7V6 (C «)2 ) _ 

oa;2 0x2 2 

( ^2 £ -l/2 + U,/6 (< ggMf) 1 s > _ e 7,/ia^). 

2 9y 2 
Next, elimination of F/2 yields 



(10) 



9x 2 l + Cfe 7 '?/ 12 (l--Cf£ 7 ''/ 12 ) 
_ £ 3/2+„ < ^i,W(^)| s > (_^-) 2 + 0( e 3 / 2 + 19 "/ 12 ). (11) 

£ 0^2 

TTie above formula results in Saffman' version of the law by Beavers and Joseph, if only the first term 
at the right hand-side is taken into consideration. For small r\, we obtain a significant deviation of 
the law by Beavers and Joseph from J 18}/ and J%jj. We are not aware of any rigorous derivation of 
a nonlinear interface law for the unconfined fluid flow coupled to the porous media flow. 
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3 Rigorous justification of the nonlinear slip law, generaliz- 
ing the law by Beavers and Joseph 

In this section we extend the justification of the law of Beavers and Joseph from [TT] to the case of 
nonlinear laminar flows. In the proofs we apply the following variant of Poincare's inequality: 

Lemma 4. (see e.g. JMS) Let ip £ V(fl E 2 ) = {ip E H 1 ^) \<p = on S e } and ip e H 1 (0 e /) such 
that ip\i X2 = e s\ = 0- Then, it holds 

\\<p\\»W<Ce 1 / 2 \\V x <p\\ ll , m , i (12) 

IMUa(nS) < Ce\\W x ip\\ L 2 {n e )2 , (13) 

L2(s) < Ce 5 / 2 \\VM LH0V ^ (14) 

<5|| 



I £8(05.*) < Ce HVxVll L * { &i-*y ( 15 ) 
3.1 The impermeable interface approximation 

Intuitive 
problem 



Intuitively, the main flow is in the fracture 0\' d . Following the approach from |11| we study the 



- e TAv° + (v°V)v° + Vp° = Fe 1 in 0f, (16) 
div v° = in £ /, (17) 



v 



°-0 on d@\ S \ ( {xi = 0} U {xi = 1} ) , (18) 



{v,p°} is 1- periodic in x u / p° dx = 0. (19) 

Therefore, as in |llj and |13) . for the lowest order approximation {v°,p } we impose on the interface 
the no-slip condition 

v° = on S. (20) 



Such choice leads to a cut-off of the shear and it introduces an error 

A solution of pre 
at E. It is given by 



A solution of problem qi6 [) -fll9 p is the classic Poiseuille flow in e { S , satisfying the no-slip condition 



v° = -e 2i -T^^(^-l)e 1 for < x 2 < e s ; p° = for < x x < 1. (21) 

2 £ £ 



Concerning the normal derivative of the tangential velocity on S, we obtain 

dv°, s _ ,F 2x 2 dv° _ ,F 



We extend v° to 02 by setting v° = for — 1 < x 2 < 0. p° is extended by to 2 . The question is 
in which sense this solution approximates the solution {v 6 ,p e } of the original problem ©-(H]). 
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A direct consequence of the weak formulation © is that the difference v e — v° satisfies the 
following variational equation 

/ e 7 V(v £ - v°)V^ dx + f ( v o d( ve ~ v °) + (^|_„0)|!L + (( v g_ v fl)v)(v e -v )V<fa! 

f div <p = / e 7 ^^ dS, V95 € W e . (23) 
Jn* Vs 9x 2 

It leads to the following result, which is a generalization of the result proved in [TTj : 

Proposition 5. Lei ms assume that (H1)-(I12) are satisfied. Let {v e ,p e } be a solution of 
and {v°,p } defined by 1121]) . Then, it holds for e < eg 

V^||V(v s - v )|| L2(Oe)i + -L||v e |U w + ||vl L2(s) + e V2-«||v« - v°|| i2(0 ^ )2 < (24) 
Proof. We test (f2"3")) with ^ = v e — v° and obtain 

l £7 | V(v e _ v } |2 dx= _ f {y e v l ){v e _ 0) M dx+ f ^M^ - „1) d 5. (25) 

Applying Lemma 2] and formula (|22p yield 

2 j«i - «o)(«I " « 2 °)^ ^1 < Ce^\\V{^ v°)|| 2 2(0 ^ )4 , 

2 (0|) 4 

Using hypothesis (HI) and above estimates lead to 

f e 7 |V(v £ - v°)| 2 dx < Ce 5+1 / 2 \\\7(^ v°)|| iW . 

We apply once more Lemma |4] and (|24|) follows. □ 

This provides the uniform a priori estimates for {v £ ,p £ }. Moreover, we have found that the 
viscous flow in E { S corresponding to an impermeable wall is an 0(e 2<5 ~ 7+1 / 2 ) L 2 -approximation for 
v e . The slip law, generalizing Beavers and Joseph's law, should correspond to the next order velocity 
correction. Since the Darcy velocity is of order 0(e 5_7+3 / 2 ), we may justify Saffman's observation 
that the bulk filtration effects are negligible at this stage. 

3.2 Justification of the nonlinear slip law 

At the interface £ the approximation from Subsection 13.11 leads to the shear stress jump equal to 
dv° F 

£ 7 — -|e = — e s . The shear stress jump requires construction of the corresponding boundary layer. 
0x2 2 



J £ 7 |^K - vo) dS\ < Ce 5+1 /2|| V(v e _ v o )!|l2( 
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X 

The natural stretching variable is given by the geometry and reads y = — . The correction 

e 

{w,p w } is given by 

-Ei- 2 A y w + e- 1 (wVy)w + e- 1 V y p w = in fi^/e U fi|/e, (26) 
div y w = in Clx/s U E/e U f2|/e, (27) 

[w](-,0) = 0; W(-,0) = and [ - e^ 1 ^-] (; 0) = 6 tM| e = l £ s on E/e, (28) 

V y w € L 2 (fi e /e) 4 and {w,^} is l/e — periodic in y±. (29) 

It is natural to rescale w and p w by setting 

w = -£ m - 7 |/%) and Pw = - £ ^(y)^. 

dv° 

Using periodicity of the geometry and independence of tt^-Ie of y, we obtain 

ox 2 

-A y /3 + V y ir= |e^ 27+2 (/3V a )/3 in lif/eU^/e, (30) 

div y /3 = in Sli/eUS/eUSl^, (31) 

[/%0)=0; H(-,0) = and [|^](.,0) = 1 on E/e, (32) 

Vy/3 e L 2 (fi7e) 4 and tt} is l/e - periodic in j/j,. (33) 

We do not use directly the nonlinear boundary layer problem (f3T)]) - (|33p . Since by (H2) we have 

F F 

5 - 2 7 + 2 > 0, we approximate {/3, n} with {/3° + — e <5 ~ 27+2 /3 1 , tt° + — e^^^n 1 }, where the new 

2 2 

functions arc given through the following problems 

-A y /3° + V.y7r° = in Qf'Ve U Q|/e, (34) 
div a /3° = in Ox/e U E/e U f2|/e, (35) 

[/3°](-,0) = 0; [tt o ](.,0) = and [^I](.,0) = 1 on S/e, (36) 

Vy^° e i 2 (^7e) 4 and {^°, tt } is l/e - periodic in yi (37) 



and 

-A^ 1 + VyTr 1 = (P Q V y )/3° in Oi' 5 /eUE/eUf2|/£, (38) 

div y p 1 =0 in Qx/e U E/e U fi|/e, (39) 

Vy/3 1 S i 2 (fi £ /e) 4 and tt 1 } is l/e - periodic in yi. (40) 

Because of the 1-periodicity of the geometry with respect to yi, problem ([34 |) -(|37 | is handled using 
Navier's boundary layer introduced in |10j . 

It reads as follows: We introduce the interface S = (0, 1) x {0}, the semi-infinite slab Z + = 
(0, 1) x (0, +oo) and the semi-infinite porous slab Z~ = U ! ^L 1 (Yp — {0, fc}). The flow region is then 
Z BL = Z+U su z-. 
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Then the following problem is considered: Find {f3 bl ,ui bl } with square-integrable gradients sat- 
isfying 

-A y p bl + V y uj u = inZ + \JZ- (41) 
div y /3 bi = in Z + U Z~ (42) 
[tf*] 8 (.,0) = and [{V/-/l}e 2 ] s (.,0)=e 1 onS (43) 
(3 bl = on U^ =1 (dY s - {0, k}), {/3 bl , w bz } is 1 - periodic in Vl (44) 

By Lax-Milgram's lemma, there is a unique f3 bl £ L 2 loc (Z B l) 2 , G X 2 (Z B l) 4 satisfying (|4"Tj) - 

(|4"4")) and w b ' € L 2 oc (iT + U unique up to a constant and satisfying (|4"Tj) . 

After [TU], [H] and [12] . we know that system (|4"Tj) - (|4"4")) describes a boundary layer, i.e. that (3 bl 
and u> bl stabilize exponentially towards constants, when 1 2/2 1 — > 00. 

Since we are studying an incompressible flow, it is useful to recall properties of the conserved 
averages. 

Proposition 6. ffTVf ). Let 

Cf= / $?(Vi,0)dyi=- [ \Vp bl {y)\ 2 dy. (45) 

JO JZ BL 

Then for every y 2 > and Vl € (0, 1), \P b \yi,y 2 ) - (Cf ,0)| < Ce~ 5y2 , for all S < 2tt. (46) 
Corollary 7. ffWj). Let 

C b J= f u bl {y u Q)d yi . (47) 
Jo 

Then for every y 2 > and y x e (0, 1), we /iave | w w (yi, y 2 ) - |< e" 27ry2 . (48) 

Proposition 8. (flOf ). Let j3 bl and uj m be defined by \4 i44\)- Then there exist positive constants 
C and 70 , such that 

\V(i bl {y u y 2 )\ + \Vu bl { yi ,y 2 )\ <Ce~^, for every ( yi ,y 2 ) e Z~ . (49) 

/3 b '(— ) is extended by zero to 2 \ fl £ . Let H be Heaviside's function. Then for every q > 1 we 
have 

||/^< £ -£(^ 0)tf (x 2 )|| Lg(02U ^ = Ce 1 ^. (50) 

Hence, our correction is not concentrated around the interface and there are some nonzero stabi- 
lization constants. We will see that these constants are closely linked with our effective interface 
law. 

As in TO] stabilization of /3°' £ towards a nonzero constant velocity C^e 1 , at the upper boundary, 
generates a counterflow. It is given by the two dimensional Couette flow d = C bl — e 1 . 

Now, after [Til], we expected that the approximation for the velocity reads 
v(e) = v° - ^e 5+1 -^ bl {- £ ) + ^s 5+1 -~'d = 

_ e 2*-7^4 ,** _ 1)e l _ ^ e 5+l-7^i(£) + ! £ «+H C fl4 e l (51) 

2 e s e s 2 y e 2 1 e 5 



9 



Concerning the pressure, there are additional complications due to the stabilization of the bound- 
ed „ $ 

ary layer pressure to C^f , when y 2 — > +00. Consequently, uj bLe — H(x2)C^ — — |s is small in fi^' 

0x2 

and we should take into account the pressure stabilization effect. 

At the flat interface E, the normal component of the normal stress reduces to the pressure field. 
Subtraction of the stabilization pressure constant at infinity leads to the pressure jump on E and 
the pressure approximation is 

p(e) = ~ e 4 (« w (-)-C£). (52) 

For the rigorous justification of the pressure approximation, leading to the pressure jump law, we 
refer to |16j . Numerical experiments, justifying independently the pressure jump are in [B]. 

We now make the velocity calculations rigorous. Let us define the errors in velocity and in the 
pressure: 

U e {x) = v e - v(e), V s {x) = f - P (e). (53) 

Remark 9. Rigorous argument, showing that U e is of order 0(e 2_7 ) 7 allows justifying Saffman's 
modification of the Beavers and Joseph law (see JE/ and JJ3/)-' On the interface E we obtain 

dvije) , _ s ^F2x 2 s F d[3\ l , r i-^ F r u . _ M, , r^J-i F 

dx2 Is - ~g 2^~ s — 2~dy~2~^' v=x/E 2 1 ~~e~ ~ 1 ^ ' ' 2' 

After averaging over E with respect to y\, we obtain the Saffman version of the law by Beavers and 
Joseph 

u f f = -eC b 1 l ^^ + 0(e 2 -^) on E, (54) 
0x2 

where is the average of V\(e) over the characteristic pore opening at the naturally permeable 

wall. The higher order terms are neglected. Nevertheless, for 7 close to 1 the Beavers and Joseph 
slip law isn't satisfactory any more. 

Next, the variational equation for {U £ ,V £ } reads 

J e^W :WLpdx + J [{U e \J)U e + (W V)v(e) + (v(f?)V)W E J ip dx 

- f W div tp dx = - [ (v(e)V)v(e)y> dx- f eyx-C\ l dS, Vtp e W e . (55) 

Note that IA £ is divergence free and the approximation satisfies the outer boundary conditions. In 
analogy with Proposition 4, pages 1120-1121, from [TT] we have 

Theorem 10. Let us suppose the hypotheses (R\)-(H2) and let U e and V e be defined by t53\). 
Then, the following estimates hold 

e\\VU s \\ L 2 {n , y + \\W\\ L 2 m . + s^WWWl^y + e x - 5 \\W\\ L2{0l , 5)2 < Ce 5 / 2 ^ (56) 

Proof. We test (|55|) by U e . Since div U E = 0, V s is eliminated from the equality. Next, arguing as 
in the proof of Proposition [5l we see that under assumptions (H1)-(H2) the viscous terms controls 
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the inertia terms. Therefore, it remains to estimate the forcing term and the interface term, coming 
from the counterflow. We have 

(v( £ )V)v( £ ) = -f (7 - s™-^ (f| - 1) + f e°+^C?£ - f )) ^tl 



2 £ fe 2 +fiA e )e dx 2 ^ £ 2s^ £ s 1 >+ 2 £ G i £ s , j 

Since W y j3 bl decays exponentially in y2 and the functions of X2 behave as X2S~ S for small X2, we 
obtain 

< c ^-2 7+ 5/2 ||w , ||L2(oe)4 (5?) 

and the leading part in the first two terms of (v(e)V)v(e) is 

Similarly, after integration by parts in f2^' and using that f3 bl is divergence free, we obtain the 
same order of e estimate as (l57l) for 



Consequently, it results in 

| / (v( £ )V)v(e)W £ dx| < Ce 35 -^+ 3 / 2 ||V^|| i2(ne) 4 (58) 
\J eUt^C bl dS\ < Ce 3 / 2 ||VW £ || L2(ne) 4. (59) 

Applying Lemma 2] yields the estimate ([55|) . □ 

Still the shear jump at the interface dominates inertia due to the counterflow. Correcting the 
shear jump term — J etpi—C bl dS is as above. The only difference is that instead of e s we have e 

and F/2 is replaced by —FC b /2. We eliminate it by modifying slightly the velocity and pressure 
corrections: 

Corollary 11. Let assumptions (H1)-(H3) hold, andU e , V 6 be defined by $53\). Let 



U x - F ~ =W- ^Cf £ 2 -^ w (-) + L^-i{C bl f%e\ (60) 
2 e 2 e° 



Then, the following estimate holds 

e\\VU^\\ L 2 {n ,y + j|W^|| L w + e^WU^Wwy + e l - 5 \\U^\\ L2{0{ ,, )2 < Cs^ 2+3S -^. (61) 
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The new shear stress iump term generated by correction (|60[) is given by — f e ~ tpi — (C\ ) dS. 
Then, the corresponding estimate (j59")l in the proof of Theorem [TUJ takes the form 

| f e 2 - s U{^{C\ 1 ) 2 dS\ < Ce b l 2 - & \\VW\\ L , [m ,. (62) 

Due to hypothesis (H3), we have 5/2 — 5 > 3(5 — 27 + 3/2 and the new error terms are less important 
than the leading inertia terms. 

Finally, we correct the inertia term effects. We note that it is multiplied by a small parameter 
e a-27+2 yj c f ji ow ^he idea from [5] and expand the solutions to the nonlinear boundary layer 
problem (|3T))) - (f3Ul) in powers of that parameter. As already explained in the beginning of the section, 

the solutions of[5D]t-(f55]) take the form {/3° + ^e 5 - 27+2 /3 1 + . . . , 7T +-e s - 2 ' y+2 7r 1 +. . . }. Furthermore, 

the 1-periodicity of the geometry in yi-dircction allows to replace j3° by (3 bl . It is similar with /3 1 . 
We recall that the leading error term for U 1,e results from (/3 W V) /3 . We introduce the boundary 
layer problem for /3 1,&i : 

-A tfi 8 1 > M +V„7r 1 ' w = G0 w V„)j0 w in Z BU (63) 
diVy/S 1 '" = in Z BL , (64) 
V/' W a 2 (^) 4 and p 1 '" 1 eL 2 oc (Z BL ) 2 , (65) 
(3 1 ' bl = on (dY s - {0, k}), and {/3 1 ' bl , tt 1m } is 1 - periodic in y x . (66) 



The forcing term decays exponentially. Following [10] . we know that the system (|63 [) -([66 [l describes 
a boundary layer, i.e. /3 1,h ' and oj 1 ' bl stabilize exponentially towards C^e 1 and C w i , when [ 2/2 1 — > 00. 
Then, the correction reads 

U 2 * =W - ^e2-T^"(£) + | £ 2- 7(c f )2 4 e i + (67) 



re 



2 7 " v e' v 2' " _11 e- 5 

In complete analogy with Theorem QJJ we prove Theorem [5J 

To obtain estimate ((8j) from Theorem[2l it is enough to note that after (|57|) . the leading remaining 
inertia terms give a contribution bounded by 



Ce 2A-+5/2-2 7 || VW 2, £ | 



Next, using hypothesis (HI), we obtain that 5/2 — 5 < 25 — 27 + 5/2. Furthermore, the leading 
order term is the shear stress jump term 



It is estimated by (|6"2"j) . which yields 
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